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Although a number of measures for quantum non-Markovianity have been proposed recently, 
it is still an open question whether these measures directly characterize the memory effect of the 
environment, i.e., the dependence of a quantum state on its past in a time evolution. In this paper, we 
present a criterion and propose a measure for non-Markovianity with clear physical interpretations 
of the memory effect. The non-Markovianity is defined by the inequality T{t 2 , to) ^ T{t 2 , ti)T{ti,to) 
in terms of memory less dynamical map T introduced in this paper. This definition is conceptually 
distinct from that based on divisibility used by Rivas et al (Phys. Rev. Lett 105, 050403 (2010)), 
whose violation is manifested by non-complete positivity of the dynamical map. We demonstrate 
via a typical quantum process that without Markovian approximation, nonzero memory effects 
(non-Markovianity) always exist even if the non-Markovianity is zero by the other non-Marovianity 
measures. 


I. INTRODUCTION 


The Markov approximation is widely used to study 
open quantum systems, leading to a reduced dynamics 
where the future state of the system depends only on 
its present state. Mathematically, a quantum Marko¬ 
vian process can be described by a master equation in 
the Lindblad forrm or equivalently by completely positive 
divisible maps[l|,l2|. Strong-coupling, finite-size environ¬ 
ments or small time scales might lead to the failure of the 
Markov approximation. Memory effects then become im¬ 
portant, and the dynamics in this case is said to be non- 
Markovian. It was found that the non-Markovianity is 
usually associated with the occurrence of revivals, non¬ 
exponential relaxation, or negative decay rates in the 
dynamics. In recent years, the features of the non- 
Markovian quantum process attracted attention in both 
theoretical and experimental studies A number 

of quantitative measures have been proposed to qua ntify 
non-Markovianity, such as the quantum channels 17 , 
the non-monotonic behaviors of distinguishability 118 - 


entanglement 221, Fisher information. correla¬ 
tion [U, the volume of states [1^, capacity [2^, the 
breakdown of divisibility the negative fidelity dif¬ 

ference [ 2 ^, the nonzero quantum discord [ 2 ^, the neg¬ 
ative d ecay rates [s^, and the notion of non-Markovian 
degree [^ . 


Among them, the measure defined by Breuer, Laine 
and Piilo [l^ is closely related to the memory effect: 
Non-Markovianity in this measure manifests itself as a 
reverse flow of information from the environment to the 
system. This back-flow of information might be a suf¬ 
ficient condition for the memory effects. It is worth 
addressing that the essence of the non-Markovianity is 
whether the future state depends on its past state or not. 
To that extent, there are no straightforward witnesses for 
non-Markovianity to date, although a great deal of ef¬ 
fort has been made to understand the non-Markovianity. 


The reason that previous measures (TtI - I^ might not 
capture directly the memory effect is the physics be¬ 
hind the definitions. The measure in Ref. jl7j| inyesti- 
gates whether a quantum channel is consistent with a 
Markovian evolution. The measure proposed by Rivas, 
Huelga, and Plenio uses divisibility as a measure of non- 
Markovianity which characterizes the noncomplete posi¬ 
tivity of the dynamical map [ 2 ^. The noncomplete pos¬ 
itivity is improved in (T7j |. but it still could not capture 
the memory effects. The other non-Markovianity mea¬ 
sures in Refs, [isl - l^ |H-[3l| are based on features dif¬ 
ferent from those mentioned before, however these fea¬ 
tures are intrinsically related to the non-complete posi¬ 
tivity of dynamical maps in quantum evolutions. There¬ 
fore, the measures so far do not directly reflect how the 
future state of a quantum system depends on its past. 
This stimulates us to consider the following questions: 
How do we directly quantify this memory effect? Is non¬ 
exponential but monotonic relaxation a Markovian pro¬ 
cess or a non-Markovian one? What is the essential dif¬ 
ference between a time-dependent Markovian equation 
and a non-Markovian time-local equation? In this paper, 
we will try to study these questions by quantifying the 
non-Markovianity directly based on the memory effect. 

Generally speaking, the conventional dynamical map 
e{t 2 , ti) is ambiguous unless the initial time of the evolu¬ 
tion is fixed. We define T{t 2 , G) (for details, see below) as 
the dynamical map transferring state from p{ti) to p{t 2 ), 
where ti (arbitrary) is the initial time of the evolution. 
The initial time means the state of the whole system is 
in a product state at this point in time. We will show 
the difference between T{t 2 ,ti) and e{t 2 ,ti) in latter dis¬ 
cussions. By this definition, T is always completely pos¬ 
itive (CP) and trace preserving (TP) for arbitrary p{ti) 
in quantum processes. The Markovian divisibility con¬ 
dition can be expressed in terms of dynamical map T 
as T{t 2 ,to) = T{t 2 ,ti)T{ti,to) because the dynamical 
map in a Markovian process does not depend on the ini- 
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tial time of evolution. In the non-Markovian regime, the 
divisible condition in terms of T is conceptually distinct 
from the conventional one since T(t2, ti) is no longer a dy¬ 
namical map connecting T{ti,to) and T(t2,to)- Instead, 
T{t2,ti) and T{t2,to) correspond to two evolutions with 
different initial times in a quantum process. Our main 
result is that the non-Markovianity in terms of T is man¬ 
ifested by inequality T{t2,to) ^ T{t2,ti)T{ti,to) which 
has a physical explanation that the evolution depends on 
its history. This is distinct from the non-Markovianity 
measure defined via c omp letely positive divisibility dis¬ 
cussed in Refs.jH, [23, HJ for the dynamical map. 

With those notations, the non-Markovianity is then 
defined as the maximum distance between T{t2,to) and 
T{t2,ti)T{ti,to) over ti and t2- It reflects the degree of 
memory effects in a quantum process. We calculate the 
non-Markovianity of a typical quantum process where a 
qubit decays into vacuum with different environmental 
spectra. The results show that a non-zero memory effect 
(non-Markovianity) always exists in this process, even if 
it is Markovian by the previous measures. Our measure 
applies to a quantum process that can describe evolu¬ 
tions starting from an arbitrary time t/, rather than an 
evolution with a fixed initial time. 

The paper is organized as follows. In Sec.II, we review 
the dynamics of open quantum systems and the concept 
of universal dynamical map. In Sec.Ill, we define two 
types of dynamical maps, A and T, with different physi¬ 
cal meanings. The non-Markovian criterion and measure 
is introduced in Sec.IV in terms of T. In Sec.V, we cal¬ 
culated and discuss the non-Markovianity of an exactly 
solvable model. Finally, we summarize in Sec.VI. 

II. DYNAMICS OF OPEN QUANTUM 
SYSTEMS 

A natural way to model an open quantum system is 
to regard it as arising from an interaction between the 
system and the environment (denoted by S and E), 
which together form a closed quantum system. The 
reduced density matrix for the system is obtained by 
tracing out the environmental degrees of freedom, i.e., 
PS = Tt^e{pse)- The total Hamiltonian for the system 
and the environment can be written as, 

H = Hs + He + Hse, (1) 

where Hs, He and Hse represents the Hamiltonian of 
the system, the environment and the coupling, respec¬ 
tively. Consider the total density matrix at ti, psE^ti), 
which undergoes a unitary evolution. The system density 
matrix at ^2 (^2 > ^i) is given by 

Psih) = T^E[U{t 2 ,ti)psE{tl)U{t 2 ,ti)^], (2) 

where U(t2,ti) = In a general case where 

the total Hamiltonian is time-dependent, U{t2,ti) can 

be expressed as U{t 2 ,ti) = with T the 


chronological operator. Eq. can be understood as a 
map £{t2,ti) connecting psiti) and psit2), namely, 

Psih) = £{t 2 ,ti)ps{ti)- (3) 

In general, such a map depends not only on the total 
evolution operator U(t2,ti) and properties of the envi¬ 
ronment, but also on the system state because psE(ti) 
may contain correlations between the system and the en¬ 
vironment 0 . Given a fixed PEiti) of the environment 
and the correlations between S and E, not all psiti) are 
allowed due to the positivity requirement of psE(ti). For 
example, if the correlation is non-zero, psif-i) can not be 
a pure state [^. Moreover, it is well known that £(^2,^1) 
may not be CP. 

A dynamical map which is independent of the state it 
acts upon is called a universal dynamical map (UDM). 
It describes a plausible evolution for any states [3 ■ Be¬ 
ing TP and CP, a UDM can be expressed in the Kraus 
representation, 

Ps{t2) = £{t2,ti)ps{ti) 

= ^ A:„(t2,ti)ps(u)A:;f(t2,ti), (4) 

a 

where ^ai^2,ti)Ka{t2,ti) = I. It turns out that a 
dynamical map is a UDM if and only if it is induced from 
the initial condition psE(ti) = Ps(U) G pE(ti), where 
PEiti) is fixed (Independent of the system state) for ar¬ 
bitrary psih) p. 

Notice that in a Markovian or non-Markovian quan¬ 
tum process, an evolution could start not only at t = 0, 
but also at the other times, say, t = tj (ti > 0) where we 
refer to tj as the initial time of the evolution. Starting at 
ti implies that the total state of the system and the envi¬ 
ronment satisfies conditions at t/. In this paper, we focus 
on quantum processes where the evolutions start with an 
initial product state of the system and the environment, 
i.e., psE^ti) = ps{ti) ® Psiti) where psiti) is fixed and 
Psiti) is arbitrary. In other words, we study models de¬ 
scribing a system in arbitrary states interacting with the 
environment in a fixed initial state. This consideration 
is reasonable as it means that the system and the envi¬ 
ronment are independent at the beginning ti and then 
start their evolution. Typically, psiti) = ps is station¬ 
ary under a time-independent environment Hamiltonian 
He, i.e., pssiti) = psiti) ® Pe for any ti. For example, 
the equilibrium state or the vacuum state of the environ¬ 
ment is often chosen as an initial environment state to 
study the reduced system dynamics. For a general (non¬ 
stationary) initial environment state, it can be given by 

psiti) = ( 5 ) 

where we have assumed that before an evolution starts, 
initial state of the environment undergoes free evolution 
driven by He [Hsit)]. 

With the factorized initial state, the dynamical maps 
starting from any time ti are UDMs characterized by the 
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fixed initial environment state and the global unitary evo¬ 
lution (total Hamiltonian), rather than the initial system 
state. Consequently, the quantum processes are supposed 
to have fixed non-Markovianity as a function of pE{ti) 
and H(t), regardless of the system states. On the other 
hand, the quantum process is universal because it can 
describe the evolution of any initial system state. In the 
case of psiti) = Pe and H{t) = H, which is the usual 
assumption used to investigate the dynamics of an open 
quantum system, we expect a time-independent value of 
non-Markovianity of the quantum process. Otherwise, 
the non-Markovianity may depend on the time interval 
that we are interested in since the quantum process is 
not time- homogenous. 


III. DYNAMICAL MAPS WITH AND 
WITHOUT MEMORY 

The reduced dynamics of an open quantum system 
can be represented by dynamical maps. In particular, 
quantum non-Markovian behaviors are often discussed 
with the help of the dynamical map e{t2,ti) that trans¬ 
forms the state from ps(ti) to ps{t2)- However, in general 
the physical meaning and properties of are not 

clear unless we know the total density matrix psE{ti), 
or, from another point of view, the initial time of the 
evolution. For instance, let to < ti < t2', when the evo¬ 
lution starts at U, which means psE{ti) = Ps(ii) ® Pe, 
the map e(t2,ti) is a UDM. In contrast, when the evolu¬ 
tion starts at a fixed time to, the total density matrix at 
h is psE{ti) = U{ti,to)ps{to) 0and then 
£(^2,^1) can be understood as an intermediate map that 
is in general not a UDM. Obviously, e(t2, U) has different 
meanings and properties for different initial times of evo¬ 
lution. This is crucial for understanding non-Markovian 
characters since it tells us that the evolution of the system 
after h [determined by e{t2,ti)] is related to its history 
(from to to ti). To clarify this point, we introduce two 
types of dynamical maps denoted by A and T, respec¬ 
tively. We will describe the two dynamical maps in the 
following. 


A. Type I: dynamical map with memory 

is defined as the quantum dynamical map 
transferring the state from ps(ti) to ps(t2) when the evo¬ 
lution starts at a prescribed time to (to < ti < ^2), i.e., 
t/ = to. When ti > to, A(t2,ti) can be understood as an 
continuous map. The physical meaning of A(t2,ti) can 
be expressed as 

Psih) = A(t2,ti)ps(ti) 

= Tr£;[C/(t2,ti)ps£;(ti)C/(t2,ti)i], (6) 

where psEih) = t/(ti,to)ps(to) 0 P£(to)Ct(ti,to)'i and 
Ps(to) = Ps{ii) defined by Eq.®. psE{ti) is in general 


not a product state with fixed environment state due to 
the system-environment interaction, this means A(t2,ti) 
is not a UDM in general. The history of the system 
from to to ti is encoded in psE{ti)- We will refer to 
this type of map as dynamical map with memory. This 
interpretation of the dynamical map is frequently used to 
investigate the non-Markovinity of quantum dynamics. 

In particular, when the starting time of the map is the 
initial time of the evolution, i.e., U = to, A(t2,to) is a 
UDM represented by, 

Ps{t2) = A(t2, to)ps(^o) 

= "^T^E\U{t2,to)ps{to) ® PE{to)U{t2,to)^\ (7) 

which can map any system state to another physical 
state. 


B. Type II: memoryless dynamical map 

Next, we define T(t2,ti) as the quantum dynamical 
map transferring the quantum system from state p(ti) 
to p{t2), where ti is the initial time of the evolution, i.e., 
tj = ti. The initial total state is psE{ti) = Ps(U) 0 
PE{ti) with fixed pE{ti) defined by Eq.® and arbitrary 
ps{ti). Then, T{t2,ti) can be expressed as, 

Ps{t 2 ) = T{t 2 ,ti)ps{ti) 

= TrE[U{t2,ti)ps{ti) 0 pE{ti)U{t 2 ,ti)'^]. (8) 

T{t2,ti) is always a UDM in contrast to A(t2,ti). As 
PE{ti) = PE{ti) is independent of the system’s history, 
the map T is memory less. We remind the reader that 
in a quantum process, T{t2,ti) and T{t2,t[) {ti ^ 
correspond to different evolutions. 

In particular, when the total Hamiltonian H is time- 
independent and the initial environmental state is sta¬ 
tionary, we have 

Ps{t2) = T{t2,ti)ps{ti) 

= Tr£;[e-A^(‘^-*i)ps(U)0pseA"(*=-‘^)].(9) 

One can easily observe that the map T{t2,ti) is time- 
homogeneous, i.e., T{t2,ti) = T{t2 — ti, 0 ). This feature 
was discussed by Chruscihski and Kossakowski in Ref. 
Q , leading to the existence of an initial time tj in a non- 
Markovian time-local master equation, 

PS = 7 ^{t - ti)ps, ( 10 ) 

such that T{t2,ti) = C.(r-t,)dr ^ C(r)dr ^ 

T{t2 — ti, 0 ). Eg. lflUl) provide a full description of the 
quantum non-Markovian process where evolutions start¬ 
ing at any time t/ are included. The initial time tj is usu¬ 
ally ignored in literature provided that one only considers 
evolutions starting from t = 0 , leading to p = C(t)p which 
may have the same form of a time-dependent Markovian 
equation with non-negative decay rates. However, the 
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existence of the initial time tj is the essential feature of 
a non-Markovian time-local equation compared with a 
time-dependent Markovian equation. It provides an evi¬ 
dent signature of the memory effect since the dynamical 

map e{t2,ti) = Te-^'i sending psiti) to ps{t2) 

relies on the initial time ti (tj < ti < <2), he., the his¬ 
tory of the system, while in a time-dependent Markovian 
equation, e(t2,ti) = Te-^o is uniquely defined. 

When <1 is the initial time of the evolution to, we have 
A(t2,^o) = T{t2,to) by definition. The key difference 
between the two types of dynamical maps A(t2,ti) and 
T{t2,ti) exists when ti is not the initial time. When 
ti to, in general A(t2,ti) ^ T(t2,ti) since psEiti) is 
in general not ps(ti) ® psiti) with fixed pEiti). Given a 
quantum process characterized by the total Hamiltonian 
and the initial environment state, the map T is well- 
defined and straightforward to calculate compared with 
A. With these notations, the non-Markovianity can be 
defined by the use of the memoryless dynamical map in 
the following section. 


IV. NON-MARKOVIANITY BASED ON 
MEMORY EFFECTS 

In this section, we give a condition as well as a mea¬ 
sure for quantum non-Markovinity quantifying the degree 
of the memory effect. Let us start from the quantum 
Markov process. A quantum process is called Markovian 
if the corresponding dynamical map, which we denote by 
Am in this paper, is divisible, i.e., 

kM{t + t) = KM{T)AM{t) ( 11 ) 

where all maps are UDMs (CP and TP) defining a 
one-parameter semigroups. Equivalently, the quantum 
Markovian process can be described by a master equa¬ 
tion in Lindblad form, [s^ 

PS = 

= -i[H,Ps\ + y^^la{VaPsVj - \VaVaPS 

OL 

-\psVjV^), (12) 

with 7 q > 0 . The evolution of the system is given 
by psit2) = AMit2 - ti)ps{ti) where Auit) = e^‘. 
More generally, when the generator C varies with time 
due to the change of external conditions, the process is 
called time-dependent Markovian which satisfies a time- 
dependent master equation in Lindblad form, 

PS = C{t)ps 

= -i[Hit),ps] +'^'rait)[Va{t)psVa\t) - 

a. 

^Vj{t)Vait)ps - ^psVJit)Va{t)] ( 13 ) 


with 7a (t) > 0 . The dynamical map becomes 

AM{t2,ti) = which is inhomogeneous but 

still a UDM, and the evolution of the system can be de¬ 
scribed by psit2) = AM{t2,ti)ps{ti). The divisibility 
condition can then be written as 

AM{t 2 ,to) = AM{t 2 ,tl)AM{tl,to) ( 14 ) 

where each dynamical map is a UDM. We will discuss 
in the following the time-dependent Markovian process 
without loss of generality. 

In a quantum (time-dependent) Markovian process, 
the dynamical map transferring ps(ti) to ps{t2) is 
uniquely given by the UDM AM(t2,ti) regardless of the 
system’s history before ti. We remark that the map 
AM{t2,ti) is a UDM if and only if it is induced from 
a product state, i.e. 

Ps{t2) = AM{t2,ti)ps{ti) 

= TTE[U{t2,ti)ps{ti) ® PE{ti)U{t 2 ,ti)^]il 5 ) 

where pEi{ti) is a fixed state and psiti) is arbitrary. 
The evolution of the system at any time U in a (time- 
dependent) Markovian process can be understood by ex¬ 
amining Eq. ITSl) . It is interpreted as the collisional model 
as reviewed in Ref.js^. Particularly, when the total 
Hamiltonian and the initial environmental state are time- 
independent, the map AM(t2,fi) can be understood as 

Psit2) = AMit2,ti)psiti) 

= AM(t 2 — U)ps(fl) 

= psiti) 0 (ig) 

which depends only on the difference t2—ti, correspond¬ 
ing to a homogeneous Markovian process described by 
Eg.lfTTl) or Eg.dT^. 

In a real Markovian quantum process, the total state 
of the system and the environment may not remain per¬ 
fectly factorized as in Eq. m, and eit2,ti) is not a UDM 
in general. From this point of view, the exact dynamics 
of an open quantum system is never perfectly Markovian 
0 . However, if psEit) ~ psit) 0 PEit) with fixed p£;(t) 
(independent of the system state) where the correlation 
between the system and environment does not affect so 
much the system’s dynamics, a Markovian model is still 
a good approximate description. 

An important observation is that the evolution from 
any time ti in a quantum (time-dependent) Markovian 
process can be interpreted as Eq. (HH) where the total 
state psEiti) remains a product state and the environ¬ 
mental state is fixed for arbitrary psiti). No matter the 
initial time of an evolution is ti or to (or any time be¬ 
fore ti), the dynamical map transferring psiti) to psit2) 
is always given by AM(t2,fi). Therefore, the Markovian 
dynamical map Am can be interpreted as both A and T 
according to our definitions in Sec.HI, i.e., 

AMit 2 ,ti) = A(t2, ti) = T(t2, ti) 


( 17 ) 
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for any to < ti < t2. Consequently, in a Markovian pro¬ 
cess, the divisibility condition Eg. (ITU) can be expressed 
both in terms of A and T, respectively, 

A(t 2 , fo) = A(f 2 , ii)A(ti, to)j ( 18 ) 

T(t 2 ,to) = T{t 2 ,ti)T(ti,to)- ( 19 ) 

Any violation of the above two divisibility conditions 
will be a manifestation of non-Markovianity. Interest¬ 
ingly, the violations of Ea. dT^ and Ea. dT^ have different 
physical interpretations which can lead to different cri¬ 
teria and measures for quantum non-Markovianity. The 
violation of Eq. ( 1181 ) is usually manifested by the non- 
complete positivity of the intermediate map A{t2,ti) in 
evolutions starting at to, which is the ultimate reason 
behind the behaviors of different quantities (such as the 
trace distance, entanglement, correlation and so on) in 
Refs . [ 1814 ^ [23, [i^- The violation itself also accounts 
for the measures in Refs. [H, 113 , lU. In contrast, 

the dynamical maps in Eq. dTOl) are all UDMs, i.e., CP 
and TP, and the violation of the second divisibility con¬ 
dition Eg. ( 1191 ) is manifested by the inequality 

T{t 2 ,to) ^ T{t 2 ,ti)T{ti,to)- ( 20 ) 

A quantum process is non-Markovian if there exists t2 > 
ti > to such that the inequality Eg. ( 1 ^ 01 ) holds. This 
criterion is conceptually different from others. It has a 
clear physical interpretation in terms of memory effects, 
which we will explain in the following. 

Consider the left-hand and right-hand sides of Eo. d^ 
acting on an arbitrary state ps(to), respectively, as shown 
in FiglU On the left-hand side of Eo. d^ . ps{to) is 
mapped to ps{t2) by T{t2,to) in evolution A, while on 
the right-hand side, ps{to) is firstly mapped to ps{ti) 
by T{ti,to) in evolution B. Then, as an initial state, 
ps{ti) is mapped to ^^(O) by T(t2,ti) in evolution C, 
which starts at ti with psE{ti) = Ps{ti) ® pE(ti)- If 
T{t2,to) ^ T{t2,ti)T{ti,to), there exists PS (to) such that 
Ps{t2) ^ p'si^'i)- We remind the reader that in both 
evolutions A and C, the system state at ti is ps{ti), 
but with different histories: in evolution A, psiti) has 
a history from to to ti which is taken into account by 
PSEih) = U{ti,to)psito) ® PE{to)U{ti,to)'^. Neverthe¬ 
less, in evolution C, ps(ti) serves as an initial state 
without any history before ti. Therefore, the fact that 
Ps{t2) Ps(^2) is a direct manifestation of the mem¬ 
ory effect where the future evolution (after ti) of the 
system depends on its history (from to to ti). This 
is the fundamental property of non-Markovianity. If 
T(t2,to) = T{t2,ti)T(ti,to) provided t2 > ti > to, then 
Psit2) = p'sit2) for any ps(to), we say that the future 
state of the system depends only on its present state, 
i.e., the process is Markovian. 

From the environment side, the physics of Eg. ( 1 ^ 01 ) can 
be interpreted as follows. At the end of evolution B, 
the total state of the system and the environment is 
PSE{ti) = U{ti,to)ps{to) ® PE{to)U{ti,to)'^ . Then, at 
the beginning of evolution C, the environment is initial¬ 
ized by T{t2,ti) such that psE{ti) -)> ps{ti) ® PE{ti), 


T{t2,to) 

I t 

Ps{to) ->- Ps{ti) - Ps{t2) 


PS {to) 


T{ti,to) 


Ps{ti) 


T{t 2 , ti) 


p's{t 2 ) 


FIG. 1: Schematic illustration of the non-Markovnian inequal¬ 
ity Eq.((20l). 


where PE{ti) is the initial environmental state at ti de¬ 
fined in Eq.([ 5 ]). The terminology initialize means the 
system information acquired by the environment in the 
time interval [to, ti] is erased at time ti . Note that the ini¬ 
tialization never happens in evolution A. After ti, if the 
future system states in evolutions A and B are different, 
that means the environment ’’remembers” the history of 
the system [encoded in the psE{ti)] and the future sys¬ 
tem is affected by this kind of memory. 

We define the non-Markovianity as the maximum dis¬ 
tance between T{t2,to) and T{t2,ti)T{ti,to) over all ti 
and t2 (We may investigate a quantum process with a 
fixed to), 

Nm = iaia.xD{T{t2,to),T{t2,ti)T{ti,to)) ( 21 ) 

Here D denotes some distance measurement be¬ 
tween the two quantum dynamical maps T{t2,ti) and 
T{t 2 ,ti)T{ti, 0 ). (Here the concatenation of the two 
UDMs T{t2,ti) and T{ti,to) is also a UDM). This def¬ 
inition allows us to quantify non-Markovianity directly 
through dynamical maps without optimization of quan¬ 
tum states [l 3 |. It can be understood as the maximal 
deviation of the divisibility condition Eq. ((T^ . When 
Nm —f 0 , a quantum process loses its memory effects 
and becomes Markovian. 

To choose a distance measure between two dynamical 
maps for our non-Markovianity, we remark that a dy¬ 
namical map A is isomorphic to its Choi-Jamiolkowski 
matrix [ill, defined as pA = A 0 I(|' 0 )(' 0 |)- Here I 
is the identity map, and |' 0 ) = ^ max¬ 

imally entangled state of the system and an ancillary 
system. The Choi-Jamiolkowski matrix has been used in 
Ref.[^ to quantify the non-complete positivity and mea¬ 
sure the non-Markovianity. Meanwhile, it is often used 
for measuring the fidelity or distance between two quan¬ 
tum channels (typically a general channel and a unitary 
one) ( 33 - [^ . Here we can easily measure the distance be¬ 
tween two dynamical maps Ai and A2 through the trace 
distance of their Choi-Jamiolkowski matrices, i.e., 

D(Ai, A 2 ) =-||pAi - PA 2 II (22) 

where ||A|| = Tr(vG 4 UA) is the trace norm of an opera¬ 
tor A. The good properties of the trace distance can be 
taken advantage of for measuring the distance of dynam¬ 
ical maps [S^. Finally, from Eq. d^Tl) and Eq. (E^ . the 
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non-Markovianity is given by 

Nm = max-||pT(t 2 .io) “ PT{t 2 ,ti)T{ti,to)\\ ( 23 ) 

ti,t 2 Z 

which naturally gives a finite value of non-Markovianity 
satisfying 0 < Nm < 1 for any quantnm process without 
normalization. 

Given a theoretically described quantum process, the 
dynamical map T{t2,ti) is always CP and describes a 
physically plausible evolution. The determination of 
T{t2,ti) (or its Choi-Jamiolkowski matrix) is straightfor¬ 
ward as long as the evolution starting from ti is known, 
regardless of how the process is described. In experiment, 
T{t2,ti) (pT(t2,ti)) could be determined through quan¬ 
tum process tomography (quantum state tomography) 
[3^. When the complete information about T{t2,ti) is 
unavailable, ps{t2) ^ can be used as a sufficient 

condition of non-Markovianity which is easy to verify 
both theoretically and experimentally. Moreover, for any 
observable A, TT:[Ap{t2)] ^ Tt:[Ap'[12)] is sufficient for 
non-Markovianity. 

V. EXAMPLE 

In this section, we illustrate how our measure can be 
calculated with a typical quantum process. The model 
describes a two-level system (denoted by S) decaying into 
its environment (denoted by E), which is initially in the 
vacuum state. This model is exactly solvable and exten¬ 
sively discussed to study the non-Markovian behaviors. 
The total time-independent Hamiltonian is written as 

H = Hs A He + Hse ( 24 ) 

where Hs = Wocr+tr”, He = are the free 

Hamiltonians of S and E, and Hse = o'-!-+ 
9 k^k denotes the coupling between the qubit and 
the environmental modes. In the case that psE(fi) = 
ps( 0 ) 0 Pe where ps( 0 ) is arbitrary and pE = | 0 )£;( 0 |e is 
the vacuum state, the evolution starting at 0 in the inter¬ 
action picture) can be expressed in terms of dynamical 
map T as 

ps{t) = T(t, 0 )ps( 0 ) 

|c(t)ppsii( 0 ) c(t)psi2(0) \ 

c(t)>S 2 l( 0 ) l-|c(t)|Vsii( 0 ) j 

where the function c(t) satisfies 

c(t) = — f dt'f{t — t')c{t') ( 26 ) 

Jo 

with the correlation fnnction f(t — ti) = 

f da;J(w)ed‘^'>-'")d-L)[i| ^ 

Since the total Hamiltonian is time-independent and 
the initial environmental state pE is invariant under 
the total Hamiltonian, the map T is time homogeneous 


according to the discussion in Sec.2, i.e., T{t 2 ,ti) = 
T{t2—ti,0). Thus, the evolution starting at ti {psE{ti) = 
Ps{ti) ® Pe) can be easily obtained as 

Psih) = T{t 2 , ti)ps{ti) =T{t 2 — ti, 0)ps{ti). (27) 

Alternatively, the evoltuion starting at 0 can be described 
by the following master equation 

PS = C,{t)ps 

2 \ 

= - 2 S{t)[a+a~,ps]+-t{t){(T-psa"^ - -(T+a~ps 

-^PSCr+cr”) (28) 

where S{t) = —2Im[^] and "/(t) = —2Re[^]. Accord¬ 
ing to Ea. (l28)) . T{t,0) is given by T(t,0) = Tedo 
We stress here that Ea. d^ can only describe evolutions 
starting at 0 rather than an arbitrary time tj. Instead, 
the master equation describing the full quantum process 
contains the initial time tj Q, 

PS = - ti)ps 

= ~ ^Ps] + lit - ti)(a~ps<J~'' 

-^a+a-ps - ^pso-+cr") (29) 

such that T{t2,ti) = Te-^‘i ^ c{T)dT 

which is consistent with Ea. dTfl) . 

Assume the spectral density of the bath are Lorentzian 
and of the form Jl{oj) = where A is 

connected to the environment correlation time te by 
te = 1 /A and 70 determines the time scale of the sys¬ 
tem by Ts = 1/70- The solution of Ea. d^ is given by 
cjt) = e-^‘/2[cosh(f )-b^ sinh(f)] with d = y'A^ - 270A 
[ij. Since c{t) is real, we have s[t — tj) = 0 in Ea. d^ . 
Typically, the exact dynamics of an open quantum sys¬ 
tem is not Markovian as discussed above. In this ex¬ 
ample, when the environment correlation time te is 
much smaller than the system characteristic time rg, 
i.e., Te a Markovian model can be a good ap¬ 

proximation P, [3^ . Thus the non-Markovianity of this 
quantum process can be characterized by the parameter 
R = ^ = ^. When i? < i, the process can be described 
by Ea. (l 29 l) with — tj) > 0 . Interestingly, when we 
only consider evolutions starting at 0 , Eq. (ESI) reduces to 
Ea. (| 28 ll which has the form of a time-dependent Marko¬ 
vian equation Eq. m ( but with different meanings be¬ 
cause T{t2,ti) is undefined in Ea. (l 28 l) for H > 0 ). There¬ 
fore, the process with i? < 4 is called Markovian by pre¬ 
viously proposed measures that use evolutions starting at 
a fixed time. It seems counter-intuitive that even though 
Te is comparable to rg, for example, R = ^ = 0 . 49 , 
the process is still called Markovian. Also, other intu¬ 
itively non-Markovian models might be called Markovian 
by previous measures ( 40 j |. In contrast, we will show that 
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FIG. 2: (Color online) Evolutions of Pee with different start¬ 
ing time. The black solid line corresponds to the first evo¬ 
lution starting at 0, and the blue dashed line corresponds to 
the second evolution starting at t\. In the evolutions, the 
system has the same state at ti, however, they became differ¬ 
ent at t 2 , which is a manifestation of the memory effect. The 
phenomenon never happens in a (time-dependent) Markovian 
process. 


non-zero non-Markovianity (memory effects) always ex¬ 
ist in this exact model, even for i? < i. In addition, the 
non-Markovinianity tends to zero as i? —^ 0. Thus our 
measure reflects how valid the Born-Markov approxima¬ 
tion is in this model. 


FIG. 3: (Color online) Non-Markovianity as a function of 
log R corresponding to the Lorentzian spectral density. It is 
shown that Nm is always non-zero for i? > 0 by our measure¬ 
ment. 


obtain the Choi-Jamiolkowski matrix of T{t 2 ,to) 

PT(t2,to) = (30) 

/ |c(r2o)P 0 


0 c{t2o) 

0 0 0 0 

0 0 1-|c(t2o)P 0 

\ c(t2o)* 0 0 1 


where T 20 = t 2 — to. Similarly, the Choi-Jamiolkowski 
matrix of T(t 2 , ti)T(ti, tg) is 


We first demonstrate that when R < 0.5, the evolu¬ 
tion in the quantum process depends on its history prov¬ 
ing the process is non-Markovian. Consider that the 
first evolution starts at 0 where the system is initially 
in its excited state, i.e., ps£;(0) = |e)(e| ® pE- Then, 
the system density matrix at ti is given by ps(ti) = 
T(ti, O)ps(O). Now we assume the second evolution 
starts at ti with initial state p'g{ti) = ps(ti) such that 
P'sEi^i) — ® Pe = Ps{ti) ® Pe- At a further time 

^ 2 , we have ps{t 2 ) = T(< 2 , 0 )/ 9 s( 0 ) in the first evolution 
and p's{t 2 ) = T{t 2 ,ti)p'g{ti) in the second evolution. The 
result is visualized in Figl^ by evaluating the evolutions 
of the excited-state population pee with R = 0.4. From 
the fact that ^ 5 (^ 2 ) Psih), we conclude that the future 
states (after ti) of the system are relevant to its history 
(from to to ti) in the process. The exited-state popu¬ 
lation Pee decays monotonically and non-exponentially 
in this case. Although the decay rate 'y{t — tj) is non¬ 
negative and the revival of pee does not occur, the en¬ 
vironment is affected by the system’s history and then 
has an influence on the future evolution of the system. 
Thus the dashed line and the solid line in Fig 15] do not 
overlap. This phenomenon never happens in a (time- 
dependent) Markovian process described by Eq. (IT51) . 
where the dynamical map transferring ps(ti) to ps(t 2 ) is 

uniquely given by A-M{t 2 , ti) = Te-^'i regardless of 

the initial time. 

The non-Markovianity Nm for this quantum process 
is calculated as follows. From Eg. (1551) and Eg. (157)) . we 


PT{t2,ti)T{ti^to) 

T{t2,ti)T{ti,to) O I(|V')(V’I) 

/ |c(r2i)P|c(Tio)P 0 


(31) 


\ c(T 2 l)*c(rio)* 


0 

0 0 

0 l-|c(r2i)p 
0 0 


c(T2i)c(rio) 

0 

|c(tio)P 0 
1 


where r 2 i = t 2 — ti and rio = ti — to. According to 
Eg. (1531) . the non-Markovnianity in the case of the Lorentz 
spectrum is 

Nm = max [ ^\M{N + \/A + N'^) \ + 

‘TiQ,r2i 8 

l\M{N-V4. + N^)\ + l\MN\ ] ( 32 ) 

8 4 

where M = c(t 2 o)-c(tio)c(t 2 i), N = c(t 2 o)+c(tio)c(t 2 i) 
and T 20 = T 21 -I- Tio. Notice that this result has been 
simplified by the fact that c{t) is real. 

We calculate the non-Markovianity by numerically op¬ 
timizing the two time differences tiq and r 2 i. Nm is 
plotted as a function of log R in Figl3] where R varies 
from 10“^ to 10"*. The result demonstrates that the non- 
Markoviaity is non-zero for all i? > 0. It monotonically 
decreases with R and tends to 0 as i? — >■ 0. It is observed 
that when R < 10“^, the non-Markovianity is already 
very small {Nm < 0.006), implying that the quantum 
process is approaching Markovian and almost memory¬ 
less. Indeed, when Nm is small, the dashed and solid 
lines in Fig|5]will be very close and almost exponential. 














Then, a Markovian master equation can well describe 
evolutions starting from any time, i.e., the full quantum 
process. In this model, the non-exponential relaxation is 
a sign of non-Markovianity. 

Before closing this section, we consider the Ohmic 
spectral density with an exponential cutoff Jo(w) = 
awe~~ , where a is a dimensionless coupling strength 
and LOc is the cutoff frequency. In this case, the full dy¬ 
namics is still described by Eas. d^Tl) or ( 1 ^ . whereas the 
function c(t — ti) is complex and S{t — tj) is non-zero in 
general. For complex c(t ), the non-Markovianity has the 
following form, 

Nm = max [-\M' - y'M'^ + A{N' - 2 K)\ + 

TlO-,T21 8 

-\M' + + - 2 K)\ + -\M'\ ]( 33 ) 

8 4 

where M' = |c(t2o)P - |c(rio)c(T2i)P, N' = |c(t2o)P -b 
|c(rio)c(T2i)P, and K = Re[c(T2o)c*(rio)c*(r2i)]. 

The exact value of for the Ohmic spectrum can 

be calculated using the analytic expressions in Ref.@, 
that is, 

c{t) = -b 

^ Jo ^uj-ujo-^{ojW + JM^)^ ^ ’ 

where E(a;) = aaJc[^e““Ei(^) — I]. When the condi¬ 
tion acoc > Wo holds, Z = [1 — ^ 0 correspond¬ 

ing to a dissipationless process (c(oo) ^ 0 ) due to the 
zero spectral density for negative frequencies, otherwise, 
Z — 0 . Here w' is the solution of w' = wo — E(w') < 0 . 
The global phase in Ea. dMl) is added compared with 
the expression in Ref. @ since we are working in the inter¬ 
action picture for consistency with the Lorentzian case. 
In fact, the non-Markovianity does not depend on the 
picture we choose. The non-Markovianity for different 
a (from 10 “^ to 10 ^) is calculated with wq/wc = I, as 
shown in Fig ID Although the system dynamics for the 
Ohmic spectrum is different from that for the lorentzian 
one, especially in the strong-coupling regime, the behav¬ 
ior of Nm as a function of the coupling strength is sim¬ 
ilar. The non-Markovianity (memory effect) is non-zero 
for all a > 0 even for parameters leading to 7(t — </) > 0 
in Eg. d^ . When the coupling is weak, log(A^M) and 
log R (or log a) are in a linear relationship in both cases, 
which indicates that the quantum process asymptotically 
becomes Markovian with the decreasing of the coupling 
strength. 

VI. CONCLUSIONS AND DISCUSSIONS 

We present a criterion and propose a measure for non- 
Markovianity of quantum processes. The measure di¬ 
rectly quantifies the degree of memory effects, i.e., how 
much the future state of a system depends on its past. To 


Nm 



FIG. 4: (Color online) Non-Markovianity as a function 
of logo corresponding to the Ohmic spectral density with 
a;o/<^c = 1. Similar to the result in FIG.®, Nm is always 
non-zero for q > 0. 


construct the measure, we introduce a universal dynam¬ 
ical map (UDM) T{t2,ti) which corresponds to an evo¬ 
lution starting from ti. In contrast to A(t2,ti), T{t2,ti) 
is well defined and always CP as well as easy to calcu¬ 
late in a quantum process. The Markovian divisibility 
can be expressed in terms of T in Eq. (fT^ and the viola¬ 
tion of it is simply manifested by the inequality Eq. ( 1201 ) , 
which has a clear physical interpretation as the memory 
effects. We define non-Markovianity as the maximal vi¬ 
olation of Ea. (|T 9 l) for all times. Unlike the previous pro¬ 
posed measures which focus on evolutions starting from 
a fixed time, our measure applies to a quantum process 
where evolutions can starting at an arbitrary time. 

One important result of our work is that a quantum 
process may have a memory effect even if it is called 
Markovian by other measures. Thus, the previously pro¬ 
posed criteria is not equivalent to the memory effect. Be¬ 
sides, we demonstrate that in a non-Markovian process, 
the dynamics may be still described by a time-dependent 
master equation in Lindblad-like form with non-negative 
decay rates. However, the non-Markovian master equa¬ 
tion contains the initial time t/, which is essentially dif¬ 
ferent from a time-dependent Markovian equation Q. 
When only describing evolutions starting from 0 (t/ = 0 ), 
a non-Markovian master equation may have the same 
form of a time-dependent Markovian master equation, 
e.g., Eg. d^ with R < \- By observing evolutions start¬ 
ing from different times as in Figj^l memory effects can be 
revealed. Thus the negative decay rates in a master equa¬ 
tion are not necessary to describe memory effects (non- 
Markovianity). And non-exponential but monotonic re¬ 
laxation may occur in both non-Markovian processes and 
time-dependent Markovian processes. 

Our measure is in units of trace distance that satisfies 
0 < Nm < 1 for any quantum process without normal¬ 
ization. It is easy to calculate regardless of the descrip¬ 
tion of the quantum process. The optimization for quan¬ 
tum states or the knowledge of the environmental state 
is not required. When the full information of the dynam¬ 
ical map T is unavailable, the condition ps(t2) 
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\Vv[Ap{t 2 )\ ^ Tr[Ap'(t 2 )]] can be used as an witness of 
non-Markovianity, which is easy to be examined both 
theoretically and experimentally. 

Now we consider the renormalized spectral density 
with zero negative frequency components in the exam¬ 
ple J'i^{oj) = Jl{oj)0{uj) where 9{u)) is the step function. 
By numerical simulation, we find that the dynamics and 
the non-Markovianity for are different from those 

for The influence caused by a negative com¬ 

ponent strongly depends on —. When — are small, 
J£(a;) and JL(a;) leads to almost the same dynamics 
and non-Markovianity for both strong coupling (large R) 
and weak coupling. When ^ are large, the negative fre¬ 
quency of the Lorentzian spectrum alters the dynamics 
and the non-Markovianity significantly for all coupling 


strengths (even for weak couplings). The reason is that 
Jl{^) is a symmetric function with x = ojq the axis of 
symmetry and A the peak width. Thus, ^ determines 
the weight of the component of the negative frequency 
which in turn alters the correlation function and the dy¬ 
namics evidently, whereas the non-Markovinity for J^(w) 
is still non-zero for all i? > 0 and tends to 0 as i? —>■ 0. 
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